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Abstract 

A buckling formulation for anisotropic 
curved panels with variable curvature is pre- 
sented in this paper. The variable curvature 
panel is assumed to consists of two or more panels 
of constant but different curvatures. Bezier func- 
tions are used as Ritz functions. Displacement 
((■°), and slope (C l ) continuities between seg- 
ments are imposed by manipulation of the Bezier 
control points. A first-order shear-deformation 
theory is used in t he buckling formulation. Re- 
sults obtained from the present formulation are 
compared with t hose from finite element simula- 
tions and are found to be in good agreement. 

Nomenclature 


M()i G) 

displacement along axial and 
transverse directions 

!/’ 

displacement along radial 
direction 

0J' , 0 y 

rotations of normals to 
middle surface or curvatures 

■G ih - 

axial, circumferential, and 
normal coordinates 


Bezier polynomial of order v 

Hi 

Bezier control points 
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£, i) non-dimensional coordinates 

L lengt h of segment 

Si widt h of segment / 

R Radius of curvature 

Introduction 

Tin' use of composite materials for aircraft 
>rimarv structures can result in significant ben- 
efits on aircraft performance and structural cost. 
Such applicat ions of composite' materials art' ex- 
pected to result in a 50-40 percent weight sav- 
ings and a 10-40 percent cost reduction com- 
pared to conventional metallic structures. Shells 
with variable curvature are widely used for air- 
craft fuselage and wing components. The vari- 
able curvature configuration for these structures 
is due to aerodynamic and functional considera- 
tions. Hence the understanding of buckling be- 
havior of composite shells with variable curvature 
is of importance for aerospace structural design. 

The earliest work reported on the buckling 
analysis of shells with variable curvature was by 
Marguerre [ 1 ] in 1951 . He expanded the curva- 
ture in a Fourier series with respect to the cir- 
cumferential arc-length in order to more accu- 
rately represent the varying curvature of wing 
leading edge panels. His approach has been t he 
basis for subsequent analyses of shells with vary- 
ing curvature for both buckling ([*2]-[4]) and vi- 
bration problems ([5]-[ll]). In these References 
Marguerre 's approach has been used for elliptical 
cylinders or panels since the variation of curva- 
ture (1 / R) can be represented accurately by one 
term in the Fourier series. Reference [4] provides 



an excellent review for buckling ami vil rat ion of 
shells and panels with variable curvatuie. 

In addition to elliptical cylinders or panels, 
shells with varying curvature such as isotropic 
conical ([12, Id]) and torisplierical shells ([13]) 
have been analyzed for free vibration lrie et 
al. [12] list'd tilt' transfer matrix technique to 
obtain natural frequencies of truncated conical 
shells. According to Reference [12], the advan- 
tage of the transfer matrix technique is its sim- 
plicity compared to other methods which require 
considerable analytical effort and computational 
time. A review of the literat ure for conical shell 
sin'll buckling is provided in Reference [12]. Singh 
([13]) used a segment approach where a series of 
circular arc segments tangential to each other at 
the segment interface juncture is used to mode] 
the shell with varying curvature (e.g., conical and 
torisplierical shell) without any approximation 
in geometry. Trigonometric and quint ic Bezier 
polynomials are used to represent displacement 
lields in each arc segment. Smooth deformed 
surface' of the shell are obtained by imposing ( ,() 
and ( ,[ continuities at the juncture of two adja- 
cent shell segments by using the properties of the 
Bezier control points. Also boundary conditions 
are applied by simple manipulation of the Bezier 
control points. Bezier polynomials have been 
used in a segment approach for the free v bration 
of isotropic elliptic cylinders ([14]), composite el- 
liptic cylinders ([15, lb]), and a laminated con- 
ical shell ([lb]). A first-order shear-deformation 
theory has been used in References [15] and [lb]. 
Detailed reviews of the free vibration of shells 
with variable curvature are discussed in Refer- 
ences [Id]- [lb]. The method presented in Ref- 
erence's [ld]-[lb] can be applied to different she'll 
geometries, shells with anisotropic material prop- 
erties, and shells with arbitrary boundary condi- 
tions. In addition, the method can incorporate 
first-order she'ar-deformaiion theory or classical 
laminated plate theory and is computationally 
efficient. 

The present analysis method for buckling of 
anisotropic shells with variable curvature use's 
a segment approach where displacement fields 
within each segment are' represented by Bezier 
polynomials and a first-order shear-deformation 
theory is used. In general, segments can be' used 
in bot h axial and circumferential directions, how- 
ever the present implementation considers only 
segments in the circumferential direction. Conti- 
nuity of displacement at the' junctures of adjacenl 


segments are' imposed using C n and ( d conditions 
obtained from the' properties of t he Bezie'r cont rol 
points ([14]). The shell with variable curvature' is 
assumed to consist of t wo or more curved panels 
of constant curvature' which is representative of 
fuselage or wing structure's. 

The present paper summarizes the analy- 
sis approach. Two structure case's with curva- 
ture' are analyzed to demonst rate' the capabilities 
of the present analysis approach. Result from 
the present analysis are' compart'd with those' ob- 
tained from finite element analyses. 

Analytical Approach 


The coordinate system and the displacement 
directions for a noncircular shell is shown in Fig- 
ure 1. Any point in the wall of the shell is speci- 
fied by means of curviline'ar coordinate system ,v. 
y and : , where ,v is the axial coordinate' fixed to 
mid-surface, y is the circumferential coordinate 
which follows the median line of the transverst' 
cross section, and c is tin' radial coordinate nor- 
mal to both J‘ and The noncircular shell is as- 
sumed to consist of two or more segments in tin' 
circumferential direction each of constant radius. 
The normal and tangent of the two segments at 
a juncture are equal as shown in Figure I, where 
;7j = ri‘j and /] = / 

Bezier polynomials are used in the axial and 
circumferential directions to represent the dis- 
placement fields. The Bezier polynomial is given 
by 


/,(».") = 


(/- 1)! (»» -/+!)! 


U'-l) 


H-/ + I 


(I) 

where n denotes t he order of the polynomial and 
0 < v < 1. For a Bezier polynomial of order n , 
there are (u+ 1 ) control points. Any point on t he 
surface of the segment is given by a parametric 
funct ion of (lie form 


A Y 

PrAZ->l) = 535Z/r(€)/.(v)»r, (2) 

r= J s = l 

where the coordinates £ and if are defined as 


£ = / L 

>i = (y - m) / (to+i - w) ( ; 5 ) 


with 0 < £.?/ < 1, A and Y are the number 
of control points in the axial and circumferen- 
tial direction respectively, and q rs are the Bezier 
control points or coefficients. The displacement 
vector can be written as 
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by 
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where «n and r n are tin' axial and transverse 


isr/r: = 

r> X A' X V' X \ SIX! - 110 X MJOIN ((>) 

for f fl> and ( fl coni inuil ies, where I SIZE is t he 
matrix size, NSEC! is the number of segments, 
and MJOIN — NS IX! — 1 for a |>anel and 
MJOIN — NS NO for a shell, i.e., MJOIN is 
t in' number of junctures. 


membrane displacements, respectively, w is the 
normal displacement and o r and d tJ are the cur- 
vatures. Subscript j = 1 . 2, T ... (A'V). The 
control points for each degree of freedom can be 
list'd to impose boundary conditions on each de- 
gree of freedom on each segment.. 

Continuity of displacement functions along 
segment junctures are obtained by using the rela- 
tions between control points of the adjacent seg- 
ment based on ( ,{l and C l continuities. Figure 
2 shows two adjacent segments and the control 
points that are involved in the ( iil and ( ,l con- 
tinuities for the case of eleven control points in 
the axial direction and six control points in the 
transverse' directions, i.e.. A = 11 and Y — G. In 
the I th segment . cont rol points e// : r; and 17^5 are re- 
lated to control points 7,1 and 7,-j of the (/ + 1 ) tfi 
segment, where /. k ~ 1 La.G according to 


<luk 




- 71/ 

- ,S/ * l:>k + <l‘>i 

S/ -f Si. j-j 



where Sf and S/ + 1 arc the widt h of the I th and 
(/ + 1 )^ J seginenl , respectively. Using these' con- 
ditions the unknowns 71; and 7 *j,- are expressed 
in terms of 754- and 7,;*. which are slaved to the 
master control points 71, and 7 v,. 

Since the buckling analysis involves first- 
order shear-deformat ion, only ( ,() continuity is 
required in the' variational formulation. However 
the advantage of also imposing (0 continuity is 
not only to obtain a more' accurate analysis but 
also to reduce the size of the stiffness and geo- 
metric stiffness matrices when a larger number 
of segments are used to represent the shell that 
is being analyzed. Table 1 shows the size of the 
matrices with the number of segments for differ- 
ent conditions of continuities when A = 1 1 and 
Y = G. If the segments are joined to approximate 
a shell, the size of the matrices are less than that 
for a panel. The size of the stiffness and geo- 
metric stiffness matrices after assembly is given 


The linear stiffness matrices are derived from 
the strain energy which is given by 



when' Ajj is the extension at stiffness coefficient, 
matrix, H fJ is the coupling stiffness coefficient 
matrix. Djj is the bending stiffness coefficient 
matrix and (' ptl is the transverse shear stiffness 
coefficient, mat rix. Tin' strain vector is { c } and 

{<} ! = {([!■ "!ryy Ky Kjry Is: ly .:} 1 (< S ) 

The strain displacement relations an* 
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Hen* ( j and ( A arc ‘‘tracer'” coefficients used to 
implement different strain-displacement relations 
or shell theories. Accordingly when ( j = (\ = 
1. the first ap])roximat ion of Sanders- Koiter shell 
theory [ 17 , 18 ] is obtained and when (\ — J, UL 
= 0 . Love’s shell theory [ 19 ] including transverse 
shear deformat ions is obtained. Finally, when C\ 
— 0 and (\ = 0 , DonnelFs shell theory [20] in- 
cluding transverse shear deformation is obtained. 





The geometric stiffness matrix is derived 
from the work done, by the applied p re- 

buckling loading and is 

I f </ — ^ ^ r ( jr A ; L ± A if * if ,V /. 

+ KrylrySl. )dA (10) 

where the nonlinear strain components are 

1 

(<> ) V L — “ ( 1 f () \ r T w \r ) 

Uu)\L = + (“\y -^)‘) 

( 5 .r y ) ;V [. — ~ M()<y 0’Oiy + ~ ) ~ f’llu- W-0,.r 

+ i«v(«\,y -~) (11) 

In the present analysis the applied prehuckling 
loading is prescribed as a uniform in-plane stress 
state. The linear stiffness and geometric stiffness 
matrices are developed using analytical integra- 
tion rather than numerical integration for com- 
putat ional efficiency. Finally, an eigenvalue prob- 
lem is solved for determining the crit ical buckling 
load. 

Numerical Results 

Results are presented for a composite cylin- 
drical panel subjected axial compression and a 
composite wing leading-edge panel subjected to 
combined axial compression and shear. Sanders- 
Koiter ([17, IS]) shell theory is used in these stud- 
ies. Buckling loads from the present analysis are 
compared with those obtained from the STAGS 
([21]) finite (dement code. The STACKS finile 
(dement model consists of the 110 element and 
curved surfaces are approximated as an assem- 
bly of flat surfaces. The nominal ply mechanical 
properties for the composite material used are: 
h 1 1 — Id. < 5 Msi; ±22 ~~ 1 4)5 Msi; G 12 — G 13 —G 23 
= 0.420 Msi and v\-> — 0.250, and the laminate 
ply stacking sequence is [±45/0/90/ ± 45],, with 
equal ply thickness for different laminate thick- 
ness. 

Cylindrical Panel 

The first structure analyzed is a semi- 
circular (n = 180") cylindrical panel 22 . 0 -in. 
Jong, and with a radius of 40.0 inches, as shown 
in Figure 5. Tin' simply-support boundary condi- 
tions are also shown in Figure 5 . Til** cylindrical 
panel is modeled as five curved segments in the 


present analysis while the STAGS finite element 
modeled consists of 20 and 40 elements in the ax- 
ial and transverse direction respectively. Table 2 
shows the results for the curved panel subjected 
to axial compression load. 

The results in Table 2 suggest that for t = 
0.072 in. the present analysis result is 4.5 per- 
cent greater than the STAGS result while for / 
= 0.144 in. and 0.2 10 in. the STAGS analysis 
results is 1.1 percent below that of the present 
analysis. The difference between results is due 
to using t he 410 shell element of STAGS which 
do not include transverse shear deformation. The 
STAGS results are above that of the present anal- 
ysis for t — 0.144 in. and 0.210 in. since these 
two panels are thicker and transverse shear de- 
formation effects are significant. The buckling 
mode shape for the curved panel obtained from 
STAGS analysis results is shown in Figure 4. 

Wing Leading- Kdge Panel 

The wing leading-edge panel is shown in Fig- 
ure 5. It consists of three curved segments of 
radii 50.0 in., 0.150 in., and 50.0 in., respectively. 
The boundary conditions of the panel is shown 
in Figure 5 and correspond to classical simply 
support conditions. Each ply of the laminate is 
0.000-in. thick. Using the present analysis, the 
wing leading edge was modeled as combination 
of two segments for the 50.0-in. radius section 
and one segment for the 0.150-in. radius sec- 
tion. The STAGS finite element model consists 
of the 410 shell element and 50x50 el ements in 
each curved segment, and the formulation of the 
410 element is based on the classical laminated 
plate theory. Table 2 shows the results obtained 
from the present analysis and those from STAGS 
for some selected combined load cases and Figure 
0 shows the buckling load interaction curve be- 
tween axial compression and positive shear load- 
ing. 

The results from the present analyse*** are 
about 4.0 percent above those of STAGS ex- 
cept for t he case of negat ive shear loading when' 
the result from the present analysis is 6.4 per- 
cent above that of STAGS. Bet ter agreement can 
be obtained from the present method by using 
more control points in the axial direction. How- 
ever this will lead to more computational effort 
and the present percentage difference between 
STAGS and results from the present analyses is 
considered acceptable for a preliminary design. 
The buckling mode shape for the wing leading 
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edge obtained from STACS analysis results are 
shown in Figure 7, 8 and 5) for axial compres- 
sion. and positive and negative shear, respec- 
tively. The shell deformation is mostly in the 
50. 0-in. radius curved segments. For the case of 
positive shear, the shell deformation is mostly it] 
our of the 50.0-in. radius curved segment. 

Concluding Remarks 

A formulation has been developed for buck- 
ling <>f anisotropic laminated shells with vari- 
able curvature. The variable curvature is ap- 
proximated by two or more segment of constant 
but dilh n nl curvatures. Bezier polynomials are 
uvd ;i> Bit/ functions in the structural axial and 
1 ran>\ < lx directions. Displacement ( ( -° ) and 
slope M’ 1 * * * * * ) continuity an' imposed between seg- 
ue in >. Ib sulls obtained from the formulation 
ar< validated using finite element simulations. 
Bin klm^ loads obtained from finite element so- 
la i h >ii" an determined to be four percent lower 
than l h I IS( I > I l he present, analysis. 
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Table 1 Size of st iffness matrices for panel and 
shell for increasing number of segments. 


nsec; 

ISIZE 

(Panel) 

[o°] [e°.e'] 

ISIZE 

(shell ) 

[r n . r 1 ] 
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:c?o 

330 
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(sor, 
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1320 


Table 2 Comparison of buckling loads results 
for composite curved panel. 


Thickness 
t (in.) 

STACS 
( Ibs/in. ) 

Present 
analysis 
( Ibs/in . ) 

0.072 

371.55 

300.08 

0.144 

1481.08 

iir>9.i. r ) 

0.216 

3328.25 

3278.86 


Table 3 Comparison of buckling loads results 
for wing- leading edge panel. 


Loading 
condition 
Aj- A r[j 

STAGS 

(Ibs/in.) 

Present 

analysis 

(lbs/in.) 

1.0 

0.0 

304.67 

317.03 

1.0 

0.4 

194.48 

202.57 

1.0 

1.0 

106.68 

110.90 

0.4 

1.0 

124.39 

129.53 

0.0 

1.0 

-117.73 

-125.37 

0.0 

1.0 

139.31 

145.55 
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Figure 4: Buckling mode shape for curved composite panel. 








Figure 7: Buckling mode shape for the composite wing leading-edge panel in axial compression 

loading. 



Figure 8: 


Buckling mode shape for the composite wing leading-edge pane] in negative shear 

loading. 
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I i^ni' !> Ruckling ino(l(‘ shape for the composite wing leading-edge panel in positive shear loading. 
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